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Abstract 

Single-chain elasticity of polyethylene at 9 point up to 90% of stretching with respect to its 
contour length is computed by Monte-Carlo simulation of an atomistic model in continuous 
space. The elasticity law together with the free-energy and the internal energy variations 
with stretching are found to be very well represented by the wormlike chain model up to 
65% of the chain elongation, provided the persistence length is treated as a temperature 
dependent parameter. Beyond this value of elongation simple ideal chain models are not able 
to describe the Monte Carlo data in a thermodynamic consistent way. This study reinforces 



*Present address: Theoretical Study of Matter, Department of Physics, University of Antwerpen, Groe- 
nenborgerlaan 171, 2020 Antewerpen (Belgium) 



1 



the use of the wormhke chain model to interpret experimental data on the elasticity of 
synthetic polymers in the finite extensibility regime, provided the chain is not yet in its fully 
stretched regime. Specific solvent effects on the elasticity law and the partition between 
energetic and entropic contributions to single chain elasticity are investigated. 
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1 Introduction 



When a single linear chain of N segments is stretched at both ends by equal and opposite 
forces ±/, the average end-to-end vector < R > lies in the direction of the force and its 
magnitude R is related to / by an elasticity law or equation of state (EOS) that is usually 
casted in the dimensionless form 

w{R,N,T) (1) 



where fee is the Boltzmann constant, T the absolute temperature, i a characteristic length 
of the unstressed chain and w a dimensionless function. 

At very low force intensities the equation of state follows the simple linear Hooke's law 

i 

where Rq is the end-to-end vector (in the mean square sense) of the polymer in its unstressed 
state. For somehow larger forces {fRo/ksT > 2) the so called "Pincus" regime is entered^, 

fRo f R^^ 



(3) 



ksT \Ro 

where 7 = 1, 3/2 for 6* and good solvent conditions respectively Note that for 6 chains 
the two regimes merge in a single one. The Pincus regime is abandoned in favor of the so 
called "Finite Extensibility" regime (FE) when the applied force is further increased beyond 
fbk/ksT = 1 where b^, the Kuhn segment of the polymer, is a measure of the local rigidity 
of the chain. In this FE regime, the EOS depends on the microscopic details of the specific 
chain being stretched and, for long enough chains, the EOS can be generally written 

{^ = w{X,T) (4) 

where X = R/Lc and Lc is the chain contour length. 

Elasticity laws of real chains in this regime are usually discussed in terms of simple ideal 
models like the Freely Jointed chain (FJC), the Wormlike chain (WLC) and some ad hoc 



extensions (elastic FJC and elastic WLC) 0, |, |], 0. The great success of the WLC 
in describing the elasticity of double-stranded DNA and of other biopolymers suggests 
the existence of some form of universality in the FE regime also. In general the ideal chain 
models predict an EOS of the form (|D with a temperature independent right hand side 

w = w{X) 

w{X) (5) 



In eq. bk is usually treated as a constant, the force being thus proportional to the 

temperature at any fixed relative extension of the chain, in good agreement with the well 
known thermo-elasticity properties of rubbers. However, the variation with temperature of 
the size of the polymer coil in a melt |T0| implies a temperature dependence of bk = 
Rq/Lci meaning that the local rigidity of the chain is affected by the temperature. In the 
context of single chain elasticity, the ideal chain models EOS eq. (|) are used to interpret 
the experimental data on the X dependence of the force. The Kuhn segment bk and the 
contour length are considered as fitting parameters in the analysis, but the temperature 
dependence of bk cannot be inferred since experiments are usually performed isothermally 
at a unique temperature^, [12], |13[. When attention is payed to energetic aspects of 
the stretching, the temperature dependence of bk becomes a sensible point. If the chain 
EOS follows eq. (|^) with bk constant, thermoelasticity predicts that the internal energy of 
the chain does not change with stretching and equivalently, that the internal force has a 
pure entropic character (/ = fs). The same EOS eq. @ with a temperature dependent bk 
predicts that the internal energy change due to isothermal stretching is proportional to the 
corresponding free energy change over the whole range of stretching and the ratio /^//(t^ 1) 
is independent of stretching. Therefore, ideal models with temperature dependent Kuhn 
segment have the whole thermodynamic complexity, and are very promising candidates to 
interpret real chain behaviors. 

The temperature dependence of the Kuhn segment has been much investigated experi- 
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mentally for a large class of synthetic polymers M and more particularly for polyethylene 



(PE) (see ref |T0| and references therein). It has been obtained from the temperature depen- 
dence of the radius of gyration of the chain in melts, from intrinsic viscosity measurements of 
the polymer in various 6 point solvents or from thermoelasticity measurements on polymer 
networks. For PE, the chain dimension in the 6 state is found to increase with decreasing 
temperature. On the basis of Rotational Isomeric State (RIS) model calculations, this trend 
has been related to the increase of trans population (respect to gauche) of the dihedral angles 
of the chain backbone as temperature drops which leads to a local straightening of the 
chain reflected by an increase of the length of the Kuhn segment. 

On this basis it is tempting to test more deeply the validity of ideal models to represent 
real chain thermo-elasticity by measuring single chain elasticity law for synthetic polymers. 
Experimentally, the interpretation of such elasticity data is made difficult by the limited 
signal to noise ratio and the limited range of forces accessible which hamper an accurate 
data fitting |, 0. 

A possible alternative route to investigate single chain elasticity is provided by numerical 
simulation of theoretical models, either coarse grained or atomistic ^ 



In the present work we use the Monte-Carlo (MC) simulation approach on a well defined 
polymer atomistic model to establish with accuracy the chain elasticity law over a wide range 
of extension (0 < R/Lc < 0.9) including the finite extensibility regime and the corresponding 
change in the chain internal energy with stretching. Our aim is twofold: i) to assess the 
validity of an EOS of the type of eq. (^ definitely simpler than the general form of eq. 
(I); ii) to test the applicability of simple ideal models to represent MC data. To follow this 
program, we consider the synthetic polymer with the simplest microscopic structure, namely 
PE and we adopt an atomistic model in continuous space. To access both weak and high 
stretching regimes with good statistical accuracy, we exploit the configurational bias MC 
methodology in combination with biased fixed-/ ensemble sampling techniques. 
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The thermoelasticity of PE has been the subject of a previous Monte Carlo study for a 
lattice model |T^ but that work is focused on the study of the solvent effect more than on 
the EOS and the test of ideal chain models. 

The paper is organized as follows: in the next section, we review the theoretical framework 
of single chain thermo-elasticity and its formal statistical mechanics basis. With the aim 
of performing single chain simulations, we review the theoretical status of a single chain 
models meant to be representative of a chain embedded in a solvent (or in a bath of similar 
chains) for which intramolecular potentials have a free energy character. In section III, 
we give for completeness the adopted PE model, namely the united-atom model devised 
originally by A. Sariban et al. [^, |23|. We give in section IV some details about the PE 
simulations performed within the fixed force ensemble, the results of which can be used 
to get also averages in the fixed end-to-end vector ensemble. In section V, we analyze the 
results for PE chains, most of them being relative to 9 point (melt state) at 400K while a few 
data are presented for PE in good solvent. We discuss the free energy and internal energy 
evolutions with stretching in the FE regime, looking at the partition between the various 
energy contributions. We analyze the adequacy of the FJC and WLC models to represent 
the thermo-elastic behavior of PE. Section VI is concerned with some results and a short 
discussion on the effects of excluded volume on the elasticity law of PE. Section VII gathers 
our conclusions. 

2 Theoretical framework 
2.1 Thermodynamics 

The central quantity in single chain thermo-elasticity is the chain Helmholtz free energy 
A{R,T) the total differential being 

dA = -SdT + l-dR. (6) 
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The corresponding single chain internal energy is then given by 

where (3 — I/ZcbT. The internal force existing within a stretched chain with end-to-end 
vector R can naturally be split into an entropic (/^) and an energetic part (/^)according to 



where 



and 



-E 



R 



2.2 Statistical Mechanics 

The system under study is a single homo-polymer chain, with degree of polymerization N , 
which is surrounded by a bath of similar (melt case) or dissimilar (solvent) molecules (note 
that in the following N will always represent the number of skeletal bonds in the linear 
chain). To simplify notations, the dependence of single chain quantities will not be made 
explicit in the following formal developments. 

The link between the Helmholtz free energy A[R, T) and statistical mechanics is 

A{R,T) = -ksT log Zr{R,T) (11) 

where the partition function Zfi{R, T), appropriate to a fixed- i? ensemble (hence the R index 
to emphasize the nature of constraint) , is expressed as 

Zr{R, T)= j d£^5 (r^ -U-R) exp (-/3t/(r,; T, p)) (12) 

where F^^ is the set of cartesian coordinates {r_i)i=Q n N+1 point particles of the chain. 

In this equation, U{r_i) is the effective potential (or potential of mean force) characterizing 



the chain embedded in a bath consisting of the other chains of the meh or the solvent in an 
infinite dilution situation. The bath is characterized by a temperature T and melt or solvent 
density p and therefore the effective potential U{T_i) originates from an integration over the 
(infinitely large number of) bath degrees of freedomFg, namely 

exp(-/5t7(ri;T,p)) = j Caexp (-/?[/(£!, £2)) (13) 

where U{T_^,T_2) is the bare total potential of the chain + bath system. 
We further note the well known equivalence 

Zr{R,T) = Z{T)W\R,T) (14) 

where Z{T) = f dRZji{R,T) is the partition function of the free chain in the bath at 
temperature T and W^{R^ T)dR is the probability to find the end-to-end vector of that free 
chain in the range between R and R + dR. 

Any single chain property 0{R,T) (such as the tensile internal force or any structural 
quantity like a bond-alignment factor) is an average computed over the subset of configu- 
rations of the free chain which satisfy a fixed- value. More explicitly, the average of the 
microscopic quantity O, denoted as < O >r, reads 

0{R, T) =< 6 >R= Zr{R, T)-' J dT,5 {rj, - - R) 6(r J exp (^-0iT,; T, p,)) (15) 

The corresponding average 0{T) for the free chain is related to the fixed- i? averages through 
the relationship 

0(T) = j dRW°{R,T) 0{R,T) (16) 

We emphasize here the special case of the internal energy E{R, T) which can be expressed as 
an average quantity of the type of eg. ([151). Indeed, due to the free energy nature (temperature 
dependency) of the potential, one gets by combining basic eqs. (0),(0),(|l2]) and (|^) the 
expression 

B(&T)=<^>, (17) 



3 Model of polyethylene chain 

Our purpose is to extract the elasticity law of a PE model which is sufficiently realistic to 
be pertinent to that specific polymer but which avoids explicit consideration of the bath 
degrees of freedom. This is done by modeling the system of interest at the level of the U{r_i) 



function for which we adopt the Sariban model ||2^, In this model based on methylene 
"united" atoms, the carbon-carbon bonds and skeletal bending angles are both rigid with 
bond lengths set to dec = 1-54 angstrom and bending angle set to 7 = 109.7 degrees. These 
geometrical data are required to evaluate the contour length of the PE chain with C-C 
bonds according to 

L, = Ar4,sin(7/2) (18) 

The intra-molecular potential, written for a chain with N+1 methylene groups, is split into 
a sum of local contributions and non local interactions, 

^(£1) = + t/"""'"'^ (19) 

where 

Af-l N-4 
jjloc ^ ^ ytor^^^^ ^ ^ (I ^^^^ _ I) (20) 

i=2 i=Q 

contains single bond torsion contributions V^"^ of the "n-butane" type taken from the 



"CHARMM" force field P4| , |25| and 1 — 5 interactions coping with the pentane effect 0. The 
non-local potential part is built as a sum of effective pairwise interactions of the Lennard- 
Jones type between united atoms at least fifth neighbor along the polymer contour. 



N-5 N 
jjnonloc ^ ^ ^ ^ 

1=0 j=i+5 



(21) 



The attractive part of the pair interactions is modulated by a factor A allowing the solvent 
quality to be varied by adjusting the relative weight of repulsive and attractive interactions. 
Lennard- Jones parameters for the present PE model are e* = OASOlKJ/mol and a* = 3.74A. 
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For a given temperature T, our PE chain will be at 9 point for a particular value A = A*(T) 
defined as the A value which guarantees that 



lim I ((ijv -ro)^) |a=a*= jim | (i?o> 



leading to the Flory characteristic ratio Coo defined by 

1 



A=A* 



N 



(22) 



lim 



A=A* 



Coo(T) 



(23) 



^-oo {Nil 

At T=400K, it was found that A* = 0.505 ± 0.005 and = 7.9 ± 0.2 [||]. This value of 
A* differs slightly from the value A* = 0.53 quoted in reference |2^ which was obtained at 
the same temperature with much shorter chains and a slightly different n-butane torsional 



form. In the same work |23[, the temperature dependence of A* for the PE model at 9 point 
was estimated to be 

\''(T^ n 9 

(24) 



^ = 2| + o.0008A-' 

T T 



Performing an equilibrium sampling of the PE chain with the potential U{r_i) yields an 
energy at fixed R given by eq. (^^ where 



dl3U 



N-5 N 

V^"" + XI 

i=0 j=i+5 



a 
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d(3X* (a^ 



(25) 



in which, according to eq . (|2^) , d{(3X*)/df3 = 0.2. 

The internal energy of the PE chain can be rewritten as a sum of potential contributions, 



E{R, T) = E^^'iR, T) + E'^'PiR, T) + ^^E"**(i?, T) 



(26) 



where the three terms are the local potential contribution, due to torsion plus the 1-5 inter- 
actions, the repulsive and the attractive parts of non local contributions respectively. They 
are ensemble averages of the corresponding terms in eq.( 
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4 Computational methodology 

Single chain Monte-Carlo simulations of PE have been performed by the configurational bias 
Monte-Carlo (CBMC) method ||2^, in combination with reptation moves of a few methylene 
units. Such a numerical procedure is relatively simple as long as the chain is free of any 
constraint on the end to end vector. We explain below how the method works for a free 
chain and how it can be adapted for sampling a fixed-/ chain ensemble, i.e. when a pair 
of equal and opposite stretching forces are applied at both chain ends. Using end to end 
distance histogram analysis, these fixed stress ensemble samplings are shown to yield the 
fixed- i? ensemble data which are more appropriate to discuss the thermo-elastic properties 
of chains. 

4.1 Sampling in the fixed-f single chain ensemble 

For free chains (not subjected to external forces), a sequence of k bonds is grown at one 
end of the chain (the particular end being chosen at random) while a sequence of k bonds 
is suppressed at the other end. Individual dihedral angles (the only relevant variables) 
associated with the growth process are generated according to a distribution oc exp {—(3Vtor), 
the rest of the interactions influencing the acceptance rate through the Rosenbluth weight 
2|. 

When simulating chains with an external stretching force along the z axis (0,0,±/), 
the associated external potential energy can be taken into account in the dihedral angle 
generation by sampling a distribution oc exp {—j3{Vtor — zf)) where z is the projection of 
the new bond in the direction of the external force. Alternatively this external energy 
contribution can be taken into account in the Rosenbluth weight. We found that the latter 
method works well as long as the force is not too strong (/ < ksT/bk) but that it is 
advantageous to switch to the former method when the force becomes higher. 

For a PE chain of N bonds sampled by MC with acceptance probability Pace for the 
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individual CBMC/reptation moves of k bonds, the generation of a new independent config- 
uration typically requires [(A^/2/c)^P^^] MC attempted steps |T^. We did experiments for 
T=400K using A = A*(400) = 0.505 for N=256 and 512 which model 6 point conditions 



25| . A few runs in good solvent at the same temperature were performed with A = 0.01 for 
N=256. All results have been produced with k = 8 methylene units for the CBMC/reptation 
moves. A minimum of about 20.000 independent configurations have been generated in each 
experiment. Such a large sampling was motivated by the need of a high statistical accuracy 
in the histogram analysis to be described below. 

4.2 Fixed-R ensemble results 

Fixed-i? ensemble results can easily be obtained from a fixed-/ ensemble sampling using 
histograms based on the instantaneous R value of the sampled stressed chain. 

Let us first remind that the fixed-/ fixed-T ensemble is characterized by a partition 
function Zf defined as 



Zf{f,T) = I dr^exp (t/(r^) - (r^ - ro) • f)] 



(27) 



with the property that Zf{0,T) is equivalent to the partition function of a free chain Z{T) 
already discussed in eq. ([T^). Let < O >f denote the average of an arbitrary microscopic 
quantity O in that ensemble. 

Let us now introduce the probability {R,T)dR to find the end-to-end vector of the 
chain in the range {R, R + dR) in the stressed chain ensemble defined by eq. (|27|) . We have 
^/^j^^ ^ / dr^6 (r^ - ro - R) exp (^(r^) - (r^ - rp) ■ f)] ^^^^ 

The distribution of the end-to-end vector in the absence of force, i.e. W^(R,N) introduced 
in eq. (|l|), is related to (R, A^) by 

W%R,N) = ^exp{-pR-i)W^iR,N), (29) 

Zj 
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where the r.h.s. is globally / independent and isotropic in R while W^CR., N) depends on / 
and has a cylindrical symmetry along the force direction. By superposing the distributions 
W'^ (known up to some multiplicative constant using eq. (]29|)) obtained from the histograms 
of ^^■^(R, A^) established for a set of simulations at several force intensities, the W^{R,T) 
profile can be obtained over a large range of R values with relatively modest computational 
efforts. 

In order to estimate 0{R,T) as defined in eq. (^) from the sampling of the stressed 
chain ensemble, one exploits the relationship 

<dS{rN -ro-K) >f=W^{R,T) < 6 >R (30) 

This expression means that 0{R,T) =< O >r can be estimated within the fixed-/ en- 
semble sampling by a simple arithmetic average of the individual O values relative to all 
configurations with an end-to-end vector lying in the relevant box in the R histogram. 

In table [l|, we list the characteristics of the different fixed-/ experiments which are ex- 
ploited in this paper. In the last column, for each force value and for the indicated sampling 
quality, we mention the range of distances which has been exploited to reconstruct the un- 
perturbed chain W^{R,T) profile according to eg . (^) . 

5 Results 

5.1 Elastic properties of the realistic PE model 

Figure |I| shows the force extension law for a 6 chain of PE at 400K obtained via fixed-/ 
ensemble calculations. The data used are, after transformation to reduced units, the points 
(< R >j, /) obtained from individual simulation runs conducted for a specific combination 
of / and values (see table |lj). To emphasize the universal character (N independent) 
of the law, we display the data corresponding to two chain lengths N=256 and 512 in 
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reduced variables, namely we plot the reduced force fbh/ksT using the equilibrium estimate 
of bk = 14.9 angstrom |^ versus the reduced elongation X = R/Lc with Lc evaluated 
according to eg. (p!8D . The force extension law in the fixed- i? ensemble (not shown) was also 
estimated using eq.(P) performing the derivative of a best fit function joining the free energy 
data of In {PA{R, T))(to be discussed later) and the elasticity curve was found to be in close 



agreement with the fixed-f ensemble one, except in the R < Rq region [^, |T8[ where the two 
ensembles strongly differ. (Note that the size of this R/Lc pathological region decreases as 
N~^/'^ with chain length). 

The data obtained for = 256 and = 512 yield a unique elasticity curve w{X) which 
will be treated in the following as a 'quasi experimental' curve corresponding to the adopted 
'realistic' 9 chain PE model. We observe a linear regime up to ~ 35% stretching which is 
followed by a rapid increase of the force due to chain FE. The discussion on the validity of 
simple chain elasticity models, like those relative to FJC and WLC models, to reproduce the 
MC data is postponed to a later section. 

We first analyze how the internal energy of the chain and its various contributions evolve 
under stretching. According to eq.(p6l), the attractive part of the non local interactions must 
be weighted with the factor d{pX*)/dp which is quoted to be 0.2 by Sariban et al. ||23l but 
that we took equal to 0.215 after the analysis detailed in section ^ To stress the universality 
of the curves, we reduce all energetic contributions by a factor ksT for dimensional reasons 
and by a factor L^/bk to work with intensive quantities. 

In figure Q , we plot the deviations of various energetic quantities from their unperturbed 
value {R = Rq) in the same thermodynamic conditions versus X = R/Lc- Data are shown for 
two chain lengths A^ = 256 and A^ = 512. Throughout this paper, the change of a quantity 
Y with stretching will be denoted as AY{R) and defined as AY{R) = Y{R) - Y{Ro). 

For N=512, the various contributions to the internal energy are shown on the same figure 
and we note that those arising from the single bond rotational energy largely dominate all 
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other terms. The local 1 — 5 interactions appear to have a marginal influence and, more 
interestingly, the non-local interactions contribute in a negligible way, in agreement with 
the 9 point conditions imposed by the choice A = A* in our simulations. We come back 
on these solvent quality aspects in section |^ where we show that good solvent conditions 
(A << A*) yield instead, as expected, a /^-dependent additional contribution to the internal 
energy E{R,T) of the chain. 

5.2 Rationalization of PE thermo-elasticity data 

The elasticity law shown in fig.|l] for our realistic model of 6 chains of PE at 400K is a 
universal curve we will denote as w{X,T). In the same figure, we compare these data 
with the functional forms relative to the FJC and of the WLC models which we denote 
respectively as wfjc{X) and wwLci^)- 'Wpjci^) is simply the inverse Langevin function 
while the WLC curve wwlc{X) and its integral WwLci^), which are required for the free 
energy calculation, were evaluated by the numerical procedure, accurate to 1.5%, which is 
reported in eqs.(13)-(14) of ref.0. We want to emphasize that the curves in figure |^ are 
obtained using the values of Lc and bx (or ip) obtained in the MC model at equilibrium and 
are not fitted to MC data. We note that the function w{X,T) based on MC data and both 
simple models curves are all starting as ~ 3X at small X and diverge (in different ways) as 
X gets close to 1. 

The origin of elasticity in FJC and WLC models, and by extension to arbitrary models 
satisfying eq.(^), is often considered as purely entropic, a character which is commonly 
interpreted as the result of the reduction of the number of conformations as the chain extends. 
This is coherent with the functional form of the eq. (§) as long as the Kuhn segment bk 
(or equivalently the persistence length ip) is constant, i.e. independent of temperature and 
external force. Eq.(^ then expresses that / is proportional to T at fixed relative extension, 
which directly implies that the force is purely entropic (see eq.(P)). The absence of any 
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energetic contribution to the force is thermodynamically equivalent to the independence of 
the internal energy respect to chain stretching, as indicated by eq.(^. When analyzing 
together fig. ^ and fig. |l|, it appears that the variation of internal energy with stretching for 
PE cannot be avoided as soon as the finite extensibility regime is entered. 

In the following, we show that PE elasticity can be modeled over a wide range of stretching 
including a large part of the finite extensibility regime on the basis of eq. where bk 
is treated as a temperature dependent quantity. Such a temperature dependence of the 
persistence length (or Kuhn segment) has been invoked earlier in the case of the WLC 
model for DNA when it models a uniform elastic rod characterized by two parameters, 
namely the contour length Lc and the elastic bending modulus k' In this model, the 
persistence length, given by ip = n'/ksT, decreases with increasing temperature as a result 
of increasing Brownian undulations of the flexible rod. Here, for synthetic polymers, the 
temperature dependence of bk or Ip must be reflected on the unperturbed size of the free 
polymer Rq at 6 point, given the relationships Rq = LJjk = ^L^ip. 

The form of eq. @with a temperature dependent persistence length has some general 
implications on the elasticity law. If we denote by W{X) the integral between and X of 
w{X), we easily get the free energy and internal energy deviations from the unstressed state 

AA{R) = h^W{R/L,) (31) 

- ip dlniT) ^^^'^^^ - dlMf) ^^^ ^ 

Eq. ( ^2]) indicate that AE{R) and AA{R) are proportional. Using eqs. (||) and (^, the 
proportionality factor satisfies 



dln{ip) _ AE{R) _ fE 



(33) 



dln{T) AA{R) f 

which implies that the enthalpic/entropic partition of the total force is independent of 
stretching. 
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In figure eq.(|32D is tested for two chain lengths (A^ = 256 and = 512). We indeed 
find a proportionahty between the free energy change and the corresponding internal energy 
change over a stretching range going up to about 65%, the best value for the proportion- 
ality factor being dln{ip) / dln{T) = —0.42. Experimental values of dln{R^)/dT for PE in 
different 9 solvents in the temperature range (373K-463K) are negative with modulus in 
the range (1.0 ^ l.2)l{]-^K-^'^. At 400K, this correspond to values of dln{Rl) / dln{T) = 
dln{£p) / dln{T) in the range —0.48 — 0.40 in agreement with our result. Previous calcu- 
lations for PE at equilibrium performed at various temperatures on a very similar model 
lead to dln{Rl)/dln{T) = —0.5 ± 0.1|2^ in agreement with our finding based on stretching 
experiments. 

The results shown in figure |^ lead to several important conclusions. The elasticity law 
of the PE model at 6 state is well represented by a law of the type as in eq. with 
a temperature independent r.h.s., up to a relative extension of 65%. The variation of the 
internal energy in this range of elongations is provided by eq. (|3^). This elasticity law 
implies that, in the range of elongations below 65% which includes a large fraction of the 
finite extensibility regime, one has indeed constant ratios Je/ f = —0.42 and fs/ f = 1-42. 
Beyond this regime of extension, the more general form of the elasticity law given in eq. (§) 
with a temperature dependent Kuhn segment, must be considered. Note that simple ideal 
model like FJC and WLC, even with temperature dependent Kuhn segment or persistence 
length, cannot provide elasticity laws of this kind. 

The above picture on the thermo-elasticity of simple chain models such as PE reveals 
that, as the temperature decreases, the average size of the chain increases as a result of the 
relative increase of the population of trans conformers. In fig.§, we show the variation with 
stretching of the gauche population (expressed in percent) Ang{R) = ng{R) — Ug^Ro) with 
respect to the gauche population in the unstretched value ng{Ro) = 65%, at T=400K. This 
curve has been obtained by the histogram outlined in eq. (|30|) on the basis of the series 
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of fixed-/ simulations for = 512 listed in table 0. The stretching variable used for the 
abscissa is chosen to test the Abe-Flory prediction Ang{R) = k{{R/ RqY ~ 1)) with k=-0.656 
which is represented as a straight line in fig.^ |]I4[. We observe that Flory's theory predicts 
quite well the population shift up to i? ~ 3i?o (corresponding to X = _R/Lc ~ 0.45 for this 
chain length). The energetic part of the force clearly comes from the increased tendency for 
local bonds to go to lower torsional energy (trans state) as R increases. 

5.3 The superiority of the WLC model 

Data displayed in fig.|l| and fig. | show respectively the variation of the reduced force w{X) 
and the reduced free energy W{X) as a function of X = R/Lc for the two PE chain lengths 
investigated. In fig.|I|, we show the predicted curves wpjc{X) and wwLci.^) corresponding 
respectively to the two basic chain models while in fig. ^, only Wwlc{^) is indicated. It 
turns out that, for the present PE case, the intermediate range of stretching is much better 
represented by the WLC model as it matches 'experimental' data up to a relative stretching 
of ~ 85% for the force-extension law. At larger stretching, we note that the limiting behavior 
of the WLC model elasticity law (Lc — i?) oc /~^^^ gives a too slow asymptotic evolution 
toward R = L^. The same qualitative picture of elasticity has been observed in experiments 
on double-stranded DNA[^, We have however already noticed as the internal energy 
variation with stretching limits the validity of EOS of the type predicted by ideal models to 
the regime X G [0, 0.65], narrower than is inferred by the force-extension behavior alone. This 
consideration puts some doubts on the present interpretation of the experimental elasticity 
laws through ideal models like FJC and WLC. We suggest that the validity of such models 
must be supported by internal energy versus extension curves. 

In the interpretation of experimental data for DNA, the regime of extreme stretching 
{X > 0.9) has been explained invoking a stretching dependent contour length Lc{f) = 
L°(l + f/fo), where L° is the contour length of the unstretched chain and /o is the elastic 
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modulus of the chain[^, |29|, ||, to be used in connection with the elasticity law of the WLC 
model. The addition of the new parameter /o, and its possible temperature dependence, 
allows the extended WLC model to overcome its previous limitation of predicting a constant 
/e// ratio[Q. Although the qualitative feature of the departure of our MC data from the 
primitive WLC model prediction is very similar to that observed for double-stranded DNA, 
the elastic WLC model is of little relevance here since the contour length of our PE chain 
model, with rigid bonds and fixed bending angles, is constant. 

We could instead define a stretching dependent persistence length ip{T,X) = £p(T)(l + 
r(X, T)), where £p{T) is the persistence length of the unstretched chain and r(X, T) a 
generic function, to be used to represent the MC data in terms of the ideal WLC model 

-^ = '"<-^-^)- (i + r(A-.r)) P"' 

Note that the temperature dependence of F is necessary in order to obtain the observed 
deviation from a constant value for the ratios AE{X) / AA{X) and fs/f- Attempting this 
analysis with our data, we find that r(X, T) remains very close to zero up to X ^ 0.6, 
is slightly negative between X = 0.6 and X = 0.85 and then diverges to oo for larger X. 
Although in the intermediate regime, the temperature derivative of r(X, T) could be invoked 
to explain the thermoelasticity behavior of our PE chain, the divergence of r(X, T) at large 
X demonstrates the inadequacy of the WLC at very high stretching. 



6 A few remarks on the solvent effects on the thermo- 
elasticity properties 

In order to investigate the solvent quality effect on PE, we performed a series of experiments 
for N=256 at 400K, but using A = 0.01 << A*(400) with the aim to model good solvent 
conditions at the same temperature. Excluded volume forces here modify the elasticity law 
as we showed earlier ||18| . 
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In figure ^ we compare the R dependence of tlie non local part of the internal energy for 
good and 6 solvents. We already noticed in figure ^ that the sum of the repulsive and attrac- 
tive energy contributions arising from non-local interactions is R independent at 6 point, a 
behavior we now show to be the result of subtle cancellation effects. Figure |^ shows the devi- 
ations AE^'^P{R, T) and Ai?"**(i?, T) as defined in eq.(^) for the different solvent conditions. 
While at 9 point, both terms exhibit linear dependence on i?, the dependence is non linear in 
good solvent. These behaviors imply that the ratio [— Aii^^'^P(i?, T)/Ai?"**(i?, T)], shown in 
the inset of figure ^, is constant in the whole range of stretching in 9 solvent only. Therefore, 
this constant value 0.215 can be taken as an estimate of d{l3\*) / 8(3 which guarantees that 
the non local part of the internal energy does not depend on stretching. We note that this 
value used to produce data in figure ^ is in agreement with the prediction of Sariban et al. 

In good solvent, the contributions of repulsive and attractive non-local forces do not com- 
pensate anymore, whatever value is adopted for d{(3X)/d(3. Assuming that A is nearly tem- 
perature independent in good solvent conditions, d{l3X) /dj3 ^ X which provides a monotonous 
decrease of internal energy with stretching (not shown). This indicates that dominant repul- 
sive non local interactions disappear progressively as the size of tensile blobs characterized 
by self-avoiding walk statistics reduces with increasing force 



7 Conclusions 

On the basis of a realistic chain model of polyethylene with solvent mediated non local 
interactions, the thermo-elastic properties of an individual chain at 9 point and in good 
solvent conditions have been investigated by Monte-Carlo simulations sampling both fixed-/ 
and fixed- i? ensembles over a large stretching regime. 

The 9 state was previously obtained by adjusting an effective parameter in the attractive 
part of the non-local interactions so that the chain size grows as N^^'^ at large A^||2^, |23| , 
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We have shown in section 6 that for such 6 point conditions, the non-local interactions do 
not contribute to the average internal force over the whole stretching range. This property 
could be taken as an alternative criterion (using a single chain length) to search for 6 point 
conditions for a particular polymer+solvent model. In good solvent conditions, non local 
interactions do contribute to the internal force for all stretching values. 

We found that the thermo-elastic properties of the chain at 6 point around 400K in terms 
of X = R/Lc are quite well represented, up to X = 65%, by an EOS where the temperature 
and relative elongation dependencies are factorized, namely f{X,T) = kBTw{X)/ip(T), 
where ip is the persistence length of the chain determined from free chain simulations ||25||. 
From the variation of the internal energy and free energy differences with stretching, we 
inferred the value of the temperature derivative of ip, namely dln{ip) / dln{T) = —0.42 in 
fair agreement with experimental observations and theoretical predictions for PE chains at 
6 point. We point out here that this temperature variation implies that the internal force at 
any stretching has a non negligible enthalpic part which is negative because extension drives 
more and more dihedral angles to adopt the lowest energy trans state. 

Up to 65%, the X dependencies of the force, i.e. the function w{X) above (see eq. (|)), 
and the of internal energy difference (see eq. (|3^) ) turn out to be very well represented by the 
WLC model, using the contour length and persistence length of our atomistic PE model at 
equilibrium, without recourse to any adjustable parameter. By comparison, the FJC model 
appears to be much less representative of the PE finite extensibility effects. This result may 
be useful to interpret AFM stretching experiments for which the minimal force which can be 
detected is above the crossover force between the scaling regime and the finite extensibility 
regime. We stress again here that the validity of the WLC model, in the present context, 
does not imply that the internal force has a purely entropic character. 

Ad-hoc extensions of ideal chain models which have been proposed so far in the literature 
do not seem to apply to the highest stretching regime explored here (65% < X < 90%). The 
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general question of thermodynamic consistency of these "elastic" WLC and FJC models 
should be carefully analyzed before they could be considered something more than simple 
fitting functions. In particular, if stretching is measured at a single temperature, the elastic- 
ity law /(X, T) and the internal energy E{X, T) must be simultaneously considered because 
thermodynamic consistency (see eq.(^)) directly relates the change of the internal force with 
temperature (at constant X) to the change of internal energy with stretching (at constant 
T). 
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parameter. For this elastic WLC model with temperature dependent persistence length 
ip{T) and contour length Lc{T) we have the following EOS 

2ftp{T) 



knT 



WWLC{X{T)) 



which implies the following relation for the thermoelasticity 



AE 



dlnin dlnLr 



dlnT dlnT 



AA ^ Lr dlnLr 



ndlnT 



X wwLc{X{T)) 



more general than the simple eq.( 
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9 Figures and tables captions 




Fig. 1: Force-extension relationships obtained in the fixed-/ ensemble for N=256 (closed 
circles) and N=512 (stars). Error bars are smaller than the size of the symbols. Continuous 
and dashed lines represent the elasticity law of the WLC and the FJC models respectively. 
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Fig. 2: Reduced internal energy difference AEbk/LckBT versus the relative extension X — 
R/L^ for a PE chain at 9 point (400K) of N=256 (closed circles) and N=512 (stars). For 
N=512, we also show the various contributions to the reduced internal energy namely the 
torsional energy (upward triangles) the 1-5 local interactions (downward triangle) and the 
non local interaction (open circles). 



27 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

R/L, 



Fig. 3: Comparison between AA{R){£p/LckBT) (closed symbols) and 
AE{R){ep/L^kBT){dln{ep)/dln{T))-'^ (open symbols) for N=256 (circles) and N=512 
(squares) versus the relative stretching R/Lc- The value [d ln{£p)/d ln{T)] = —0.42 has 
been adopted to rescale the internal energy difference. The continuous curve represents 
WwLc{x) as defined in the text. 
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Fig. 4: Percentage of gauche population with respect to the unstretched state Aug — ng{x) — 
ng{0) in the fixed- ensemble as a function of {R/RqY — 1 for a N=512 chain at 9 point. 
The straight line is the Flory's prediction for small stretching. 
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Fig. 5: Attractive (squares) and repulsive (circles) contributions of the non local part of the 
energy with respect to the unstretched state in good solvent (open symbols) and 6 solvent 
(closed symbols) at 400K. Energies arc normalized to ksT. In the inset we show the ratio of 
minus the repulsive part to the attractive part which is constant for solvent but depends 
on stretching for good solvent. 
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Table 1: List of MC experiments in the fixed-/ ensemble for two different system sizes N at 
9 point. The last column gives the range of extensions statistically accessible, for the given 
imposed force / and number of statistically independent configurations of the chain. 



N 


imposed / 


# indep. configurations 


range of R/Lc exploited 


256 


0.0 


73100 


0.00 - 0.25 




0.2 


180400 


0.25 - 0.47 




0.5 


153500 


0.47 - 0.62 




1.0 


95400 


0.62 - 0.75 




2.0 


27000 


0.75 - 0.81 




3.0 


27400 


0.81 - 0.87 




4.0 


25000 


0.87-0.92 


512 


0.0 


36000 


0.00 - 0.28 




0.3 


35000 


0.28-0.51 




0.5 


23000 


0.51-0.60 




1.0 


24000 


0.60 - 0.71 




2.0 


19200 


0.71 - 0.81 
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Figure captions 

Figure 1: Force-extension relationships obtained in the fixed-/ ensemble for N=256 
(closed circles) and N=512 (stars). Error bars are smaller than the size of the symbols. 
Continuous and dashed fines represent the elasticity law of the WLC and the FJC 
models respectively. 

Figure 2: Reduced internal energy difference AEbk/LcksT versus the relative ex- 
tension X — R/Lc for a PE chain at 9 point (400K) of N=256 (closed circles) and 
N=512 (stars). For N=512, we also show the various contributions to the reduced in- 
ternal energy namely the torsional energy (upward triangles) the 1-5 local interactions 
(downward triangle) and the non local interaction (open circles). 

Figure 3: Comparison between AA{R)(£p/ LcksT) (closed symbols) and 
AE{R){ep/L^kBT){dln{ep)/dln{T))-'^ (open symbols) for N=256 (circles) and N=512 
(squares) versus the relative stretching R/L^. The value [d ln{ip)/d ln{T)] = —0.42 has 
been adopted to rescale the internal energy difference. The continuous curve represents 
WwLc{x) as defined in the text. 

Figure 4: Percentage of gauche population with respect to the unstretched state 
Aug — ng{x) —ng{0) in the fixed- ensemble as a function of {R/RqY — 1 for a N=512 
chain at 6 point. The straight line is the Flory's prediction for small stretching. 

Figure 5: Attractive (squares) and repulsive (circles) contributions of the non local 
part of the energy with respect to the unstretched state in good solvent (open symbols) 
and 9 solvent (closed symbols) at 400K. Energies are normafized to ksT. In the inset 
we show the ratio of minus the repulsive part to the attractive part which is constant 
for 9 solvent but depends on stretching for good solvent. 
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